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The impact parameter dependent gluon distribution of the proton xG{x,Q'^ ,\h±_\) is inves- 
tigated in a loop-loop correlation model that respects the S-matrix unitarity condition in 
impact parameter space. We find low-x saturation of xG{x,Q'^ ,\h±_\) as a manifestation of 
S'-matrix unitarity. The integrated gluon distribution xG(x, Q^) does not saturate because of 
the growth of the effective proton radius with decreasing x. 

1 Introduction 

The steep rise of the gluon distribution xG{x, Q^) and structure function F2{x,, Q'^) of the proton 
towards smah x = Q'^/s is one of the most exciting results of the HERA experiments § As 
F2{x,Q'^) is equivalent to the total 7*p cross section, a^J^p{s , Q'^) , this behavior reflects the rise 
of a^°^p{s,Q'^) with increasing cm. energy which becomes stronger with increasing photon 
virtuality . In hadronic interactions, the rise of the total cross sections is limited by the 
Proissart bound, which is a direct consequence of S'-matrix unitarity, 55^ = S'^ S = H, and 
allows at most a logarithmic energy dependence at asymptotic energies. Analogously, the rise 
of (T*°p(s,(5^) is expected to slow down. The microscopic picture behind this slow-down is the 
concept of gluon saturation: Since the gluon density in the proton becomes large at high energies 
^/s (small x), gluon fusion processes are expected to tame the growth of a!:^*p{s, Q^), and it is a 
key issue to determine the energy at which these processes become significant. 

In this talk, gluon saturation is considered in an effective loop-loop correlation model 
(LLCM) that respects the S'-matrix unitarity condition in impact parameter space and allows 
a unified description of pp, j*p, and 77 reactionsp Concentrating on 7^^ reactions, the impact 
parameter dependent gluon distribution of the proton xG{x, Q'^,\b±\) is computed and found to 
saturate in accordance with S-matrix unitarity. The presented results are extracted from Ref. 2 
where more details can be found. 



2 The Loop-Loop Correlation Model 



Recently, we have developed a loop-loop correlation model (LLCM) to compute high-energy 
hadron-hadron, photon-hadron, and photon-photon reactions involving real aniL virtual photons 
as wellp Based on the functional integral approach to high-energy scatteringPo'El the T-matrix 
element for elastic reactions with transverse momentum transfer q_\_ {t = —Qi), cm. energy 
squared s, and photon virtuality reads 

T^*p{s,t,Q^) = 2is[dH^e"^'^J^*p{s,\h\,Q^) (1) 



Jylp{s, \b±\,Q'^) = j dzid^rij dz2d'^r2\Tpyi{zi,fi,Q'^)\'^\'ipp{z2,f2)\'^ (l - 5'dd(6i, ^^i, n, 2:2, ^2)) (2) 
where the correlation of two light-like Wegner- Wilson loops, the loop-loop correlation function, 



SDD{b±,zi,n,Z2,f2) = {W[Ci]W[C2])^ With W[Q] = ^TTVexp 
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(3) 



describes the elastic scattering of two light-like color-dipoles (DD) with transverse size and 
orientation rj and longitudinal quark momentum fraction Zj at impact parameter 6j_, i.e., the 
loops Ci represent the trajectories of the scattering color-dipoles. For elastic scattering, the 
color-dipoles are given by the quark and antiquark in the photon and in a simplified picture by 
a quark and diquark in the proton. The and Zi distributions of these color-dipoles are given 
respectively by the perturbatively derived longitudinal photon wave function IV'^^ (zi, ri, Q^)^ 
and the simple phenomenological Gaussian wave function |V'p(-22) '^2)P- To account for the non- 
perturbative region of low in the photon wave function, quark masses mf{Q'^) are used that 
interpolate between the current quarks at large and the constituent quarks at small Q^. 

In contrast to the^wave functions, the loop-loop correlation function Sdd is universal for pp, 
7*p, and 77 reactionsp We have computed Sdd in the Berger-Nachtmann approachp in which 
the S'-matrix unitarity condition is respected as a consequence of a matrix cumulant expansion 
and the Gaussian approximation of the functional integrals. To describe the QCD interaction 
of the color-dipoles, we have used the non-perturbative stochastic vacuum modelU supplemented 
by perturbative gluon exchange. This combination allows us to describe long and short distance 
correlations in agreement with lattice computations of the gluon field strength correlator! and 
leads to the static quark-antiquark potential with color-Coulomb behavior at short distances and 
confining linear rise at long distancesB Two components are obtained of which the perturba- 
tive (P) component, (x^)'^) describes two-gluon exchange and the non-perturbative (NP) com- 
ponent, (x^^)^) the corresponding non-perturbative two-point interactiont3 Ascribing a weak 
(e^^ = 0.125) and strong (e^ = 0.73) powerlike energy dependence to the non-perturbative and 
perturbative component, respectively, 
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our final result for Sod reads 



2 '1 NP t \ \ f ^ P f \\ , ^ NP f \\ ( ^ P / 



Sdd = 3 cos(^-x''^(5)j cos(^-x^(s)j + 3 cos(^-x''^(s) j cos(^-x^(5)j • (5) 

The cosine functions ensure the unitarity condition in impact parameter space as they average 
to zero in the integration over the dipole orientations at very high energies. In fact, the higher 
order terms in the expansion of the cosine functions describe multiple gluonic interactions and are 
crucial for the saturation of the impact parameter dependent gluon distribution xG{x, Q^, l^il)- 
Moreover, they tame the growth of the cross sections at ultra-high energiesp 
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Figure 1: (a) The profile function {■K/a)J.f^p{s,\b±\,Q'^) as a function of the impact parameter at photon 
virtuality = 1 GeV^ for cm. energies from ^/s = 10 GeV to 10^ GeV. (b) The gluon distribution of the proton 
at zero impact parameter, xG{x,Q^, \b±\ = 0), as a function of x = Q^/s for — 1, 10, andlOOGeV^. 



3 Gluon Saturation 

Based on the close relationship between the longitudinal structure function Fl{x,Q'^) and the 
gluon distribution of the proton xG{x,Q'^) at small x, the impact parameter dependent gluon 
distribution xG{x, Q^, \b±\) has been related to the profile functionB J-y*p{s = Q'^/x, \b±\, Q^) 



QM 6x1) « 1.305-^ J^*p(0.417x,|6x|,Q2) . (6) 



Consequently, the shape of xG{x,Q'^ ,\b_\_\) is determined by the profile function which is 
a measure for the blackness or opacity of the interacting particles and, thus, gives an intuitive 
geometrical picture for the energy dependence of "f*j^p reactions: As shown in Fig. |l]a, the in- 
teracting particles become blacker and larger with increasing cm. energy y/s. At ultra-high 
energies, y/s ^ 10^ GeV for = 1 GeV^, the opacity saturates at the black disc limit first for 
zero impact parameter while the transverse expansion of the scattered particles continues. 

In the LLCM with a proton wave function normalized to one, the black disc limit is given 
by the normalization of the longitudinal photon wave function 

J^p{Q^)= fdzd\\i;,*^{z,r,Q^)\'' (7) 



and is reached at very high energies when the cosine functions in Sdd average to zero in the Zi 
and Tj integrations (^). 
low-x saturation value 



and Tj integrations @. Accordingly, J™p^(Q^) induces an upper bound on xG{x,Q'^, \b±\), the 



xG{x, QMSiI) < xG"^'^"(Q2) ~ 1.305 -^L 1 JYt%^{Q'') « , (8) 

n as a "-^ ots 

which is consistent with complementary investigationslHI and indicates strong color-field strengths 
G^j^ ~ Xj yfoTg as well. Since the black disc limit is a rigid unitarity bound for purely imaginary 
elastic amplitudes, i.e., obtained directly from the S'-matrix unitarity condition in impact pa- 
rameter space, we conclude that the low-x saturation of the impact parameter dependent gluon 
distribution xG{x, Q'^,\b±\) is a manifestation of S'-matrix unitarity. 

In Fig. the small-x saturation of the gluon distribution at zero impact parameter 
xG{x,Q'^,\b±\ = 0) is illustrated for = ig, andlOOGeV^ as obtained in the LLCM. The 



saturation occurs at very low values of x ^ 10"^" for ^ 1 GeV^, where the photon virtuality 
determines the saturation value (^) and the Bjorken-x at which it is reached. For larger Q^, 
the low-x saturation value is larger and is reached at smaller values of x. Moreover, the growth 
of xG{x,Q'^ ,\b±\ = 0) with decreasing x becomes stronger with increasing Q^. This results 
from the stronger energy increase of the perturbative component in the LLCM, = 0.73, that 
becomes more important with decreasing dipole size. In conclusion, our approach predicts the 
onset of the xG'(x, Q^, |6j_|)-saturation for > 1 GeV^ at x < 10"^*^ which is far below the 
x-regions accessible at HERA (x k, 10"^) and THERA (x lO"'^). 

Note that the S'-matrix unitarity condition together with relation (^) requires the saturation 
of the impact parameter dependent gluon distribution xG(x, Q'^, but not the saturation 
of the integrated gluon distribution xG{x,Q^). Indeed, approximating xG(x, Q^, in the 
saturation regime by a step-function, xG{x,Q'^ ,\b±\) xG^°-^{Q'^)Q[R{x,Q'^) — where 
R{x, Q"^) denotes the full width at half maximum of the profile function, one obtains 

Q-) . 1.305 ^ . , Q^RH^,Q-) 

which does not saturate because of the increase of the effective proton radius R{x, Q"^) with de- 
creasing X. Nevertheless, although xG(x, Q'^) does not saturate, the saturation of xG{x, Q^, 
leads to a slow-down in its growth towards small x. 

4 Conclusion 

We have computed the impact parameter dependent gluon distribution xG(x, Q^, from the 
profile function for 'jf^p reactions in a loop-loop correlation model (LLCM). The LLCM combines 
perturbative and non-perturbative QCD in agreement with lattice investigations, provides a 
unified description of pp, 7*p, and 77 reactions, and respects the S'-matrix unitarity condition 
in impact parameter space. As a manifestation of S'-matrix unitarity, we have found low-x 
saturation of xG{x, Q"^ , \b±\) for Q"^ > 1 GeV^ at x ^ 10~^'^ but only a slow-down of the integrated 
gluon distribution xG{x,Q^) since the effective proton radius grows with decreasing x. 
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